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1 Introduction 



Q . Scholte waves are acoustic waves propagating at a fluid/solid interface. They 

■^ ! are localized in the neighborhood of the phase boundary in the sense that 

c^ I they decay exponentially in both directions along the normal to the interface. 

Johnson \T\ established the exphcit secular equation for Scholte waves over 
_L ■ an orthorhombic crystal. In his case, the crystal is cut along a plane 0:2 = 

Ch ■ containing two crystallographic axes Oxi and Ox^; the wave propagates with 

2 ! speed V in the Xi direction; the solid is characterized by a mass density ps 

and relevant elastic stiffnesses Cn, Cu, C22, and Cqq, the fluid by a mass 
density pf and speed of sound c. The secular equation is 

"^ [ Z y C11C22 — C12 — 2C12CQQ — (C22 + Cqq)X + 2^^/C22CQQ{Cll — X){Cqq — X) 

^: /Cfifi - X 



. where 



(C11C22 — C12 — C22X) + X\/C22Cqq — 0, (1.1^ 



X = p,v\ Z = ^i^. (1.2) 

'i-4 



X 

^ , 

c^ , For instance, consider a frozen lake with a layer of ice assumed thick enough 

to be considered as a semi-infinite body. At 0.01°C under 1 bar, the density 
of water is [2j: Pj = 999.84 kg/m'^ and sound propagates at c = 1402.4 
m/s; the second line of Table 1 lists the elastic stiffnesses and density of 
ice [3]; according to (11. ip . Scholte waves propagate for this model at speed 
f 5 = 1237.6 m/s. Ice however has the special property of being transversally 
isotropic, which means that any plane containing the X3 axis is a symmetry 



plane and so the speed vs is the same for any orientation of the water/ice 
interface plane containing the 0:3 axis. 

The aim of this Letter to the Editor is to derive explicitly the secular 
equation for Scholte waves at the interface between a fluid and an anisotropic 
crystal cut along a plane containing the normal to a single symmetry plane, 
that is containing only one crystallographic axis. In effect, the crystal may be 
a monoclinic crystal with symmetry plane at X3 = 0, or a rhombic, tetragonal, 
or cubic crystal cut along a plane containing X3 and making an angle 6 y^ 
with the other crystallographic planes; the higher symmetry cases (^ = or 
transversally isotropic and isotropic crystals) are covered by (11.11) . For cases 
with less symmetries, one can turn to approximate solutions [1] as long as 
the anisotropy is weak . 

2 Equations of motion and boundary condi- 
tions 

Consider two half-spaces delimited by the plane X2 = 0; the upper one X2 < 
is filled with an inviscid fluid, the lower one X2 > is made of a monoclinic 
crystal with symmetry plane at X3 = whose relevant non-zero reduced com- 
pliances are s'^^^, S22, s[2, s'iq, S2q, and Sgg. At the interface, an inhomogeneous 
plane wave travels with speed v and wave number k in the xi direction, and 
decays rapidly in the X2 -> ±00 directions. 

In the solid, the corresponding equations of motion are written as a first- 
order differential system for the 4-component displacement-traction vector, 

^' = iN^, ^ikX2) = [UiikX2),U2ikX2),h2ikX2),t22ikX2)r, (2.1) 

where the functions Ui and tj2 are related to the in-plane mechanical dis- 
placements Ml, U2 and in-plane tractions au, (^ti through 

Mi(xi,a;2,a;3,t) = Ui{kx2y^^''^~'"^^ , ai2{xi,X2,X3,t) = ikti2{kx2y^^''^~''^\ 

(2.2) 
In (12.11) , the 4x4 matrix N is given by O [6] , 



N 



-re -1 nee "-26 

-r2 n2e ^ee 

X — f] —re — r2 

X -1 



(2.3) 



where X = psv'^ and 
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r,- = -■ 



'li 



'11 
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«j 
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(2.4) 



These equations also cover the case of a wave (12.21) travelhng in a crystal 
of rhombic, tetragonal, or cubic symmetry, with acoustic axes XYx^ and 
reduced compliances 5"^^, cut along the plane 0:2 = containing the X3 axis 
and making an angle 6 with the crystallographic XY plane (see Figure 1). 
In that case, the reduced compliances s^ • along the Xi axes are given in terms 
of those along the crystallographic axes XYx-^ by (see Ting [7|), 

s[^ = S[^ cos^ e + {2S[^ + S'^e) cos^ 6 sin^ 6 + S'^^ sin^ 9, 
s^2 = ^22 cos^ e + (25^2 + S'qq) cos^ 9 slu^ 9 + S[^ sin^ 9, 



s' 



12 



s' 



= S[^ cos^ 9 + {2S[^ + S'^e) cos^ 9 sin^ 9 + S'^^ 
= S'22 cos^ 9 + (25^2 + Sea) cos^ 9 sin^ ^ + S[^ 
= S[, + {S[, + S',, - 2S[, - S'J cos^ 9 sin^ 9^ 
= S',, + 4{S[, + S^2 - 2^12 - ^66) cos^ 9 sin^ ^. 



s'le = [2^22 sin^ ^ - 25^ cos^ 9 + (25^2 + ^ee) (cos^ ^ - sm^ ^)] cos ^ sin 9, 
^26 = [2'522 cos^ 9 - 2S[-^ sin^ 9 - (25^2 + 566)(cos^ 9 - sin^ e)] cos 9 sin ^. 

(2.5) 

Note that for transversally isotropic crystals, the following relationships hold, 
'S'li = 5'22, 5'gg = 2{S[^ — S[2), and the rotation does not affect the values 
of the compliances (s^^ = S'^j). This author [8] recently showed that for 
waves vanishing with increasing distance from the plane X2 = 0, the following 
fundamental relationships hold for any positive or negative integer power n 
of the matrix N, 



^(0) • IN"^(0) = 0, where I 
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(2.6) 



Because of the Cayley-Hamilton theorem, only three consecutive powers of 
N are linearly independent so that (12. 6 p reduces to only three linearly inde- 
pendent equations. 

In the fluid, the normal displacement and the normal stress component are 
connected, as recalled by Barnett et al. [9J, by the (real) normal impedance 
Z defined in (IL2D2, 

(722 = kZu2. (2.7) 

At the solid/fluid interface, the normal displacement and the normal 
stress component are continuous, and the shear stress component is zero. 
It follows from these boundary conditions and from (I2.IP 9. (12. 2p . (12. 7p . that 
the displacement-traction vector at the interface X2 = 0"*" is of the form, 

^(0+) = ?72(0)[a,l,0,-iZr, (2.8) 



wherea = f/i(0+)/f/2(0). 

Now the fundamental equations fl2.6p read 

Writing a as a = ai+ia2 and taking in turn n = 
linear system of equations follows, 



= -(iV")42-^'(iV")24. (2.9) 

— 1, 1, 2, a non- homogeneous 
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N* + Z'^N; 



24 



N^2 + Z^N2^ 
Z\N^\ 



;24 



(2.10) 



where N* denotes the adjoint of N. The unique solutions to the system are 
hk = Afc/A, where A = det A and A^ is the determinant of the matrix 
derived from A by replacing the k-th column with d. However, the b^ are 
linked by b\ + h\ = 463, which is the explicit secular equation for Scholte wave 
over a monoclinic crystal with symmetry plane at 0:3 = 0, 



A? 



A^ 
^2 



4AA. 



(2.11) 



As a check, the limit case of a solid/vacuum interface is examined. When 
the density of the fluid p/ is taken as zero, then by fll.2p 9 Z = 0, and so 
A = Ai = A3 = 0. The secular equation reduces to A2 = (written at 
Z = 0), that is the following quartic in X = pgv"^ [101 El E], 



X[r2n-n2eiX-r])] (X 

(l+r2)X-r/ 



r]){l+neeX)+rlX X[rl - n,eiX - r])] 

X X -7] 

2r6(X-r7) 

(2.12 



0. 



3 Examples 

Calculations for usual combinations of a solid and a fluid show that in gen- 
eral the speed of a Scholte wave is very close to the speed of sound in the 
fluid. Hence, consider water (pj = 1025 kg/m^, c = 1531 m/s at 25°C [TT] ) 
over gypsum (monoclinic, Ps and Cij in Table 1 [12]): the secular equation 
(12. lip yields a Scholte wave speed within the interval [1519 m/s, 1526 m/s] 
(depending on the orientation of the cut plane), which is within less than 



0.8% of the speed of sound in the water and beyond reasonable accuracy for 
measurements. 

Yet for certain choices, the Scholte wave speed moves away from the speed 
of sound in the fluid. One example is the combination ice/water presented 
in the Introduction. A second example is the combination of pure water 
(pj = 998 kg/m^, c = 1498 m/s at 25°C [11]) and Terpine Monohydride 
(orthorhombic, Ps and Cij in Table 1 |3]): at ^ = 0° and 6 = 90° (crystal cut 
along a plane containing two crystallographic axes) the wave propagates at 
1228.3 m/s and 1249.5 m/s, respectively; Figure 2(a) shows how the Scholte 
wave speed varies between these two extremes as a function of 6. Another 
way of separating distinctly the Scholte wave speed from the sound speed is to 
increase the pressure, and hence the speed of sound, in the fluid. Crowhurst 
|13j et al. recently measured the Scholte wave speed for Methanol over Ger- 
manium in a diamond anvil cell: as the pressure increases from 0.56 GPa 
to 2.2 GPa, so does the speed of sound in Methanol, from about 2500 m/s 
to 3500 m/s. In Table 1, the stiffnesses and density of Germanium (cubic) 
at 20° are recalled [5]; the density of Methanol is 791.4 kg/m^ at 20° [TT] . 
Figure 2(b) shows, in agreement with their results, the combined influence 
of orientation and speed of sound on Scholte wave propagation; each curve 
corresponds to a different speed of sound in Methanol, from c = 2000 m/s 
(bottom curve) to c = 4000 m/s (top curve) by 500 m/s increments. 
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Figure 1: Fluid/solid interface 
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Figure 2: Scholte wave speeds for (a) Water/Terpine interface and (b) 
Methanol/ Germanium interface, where the speed of sound in the fluid is 
[m/s]: 2000 (bottom curve), 2500, 3000, 3500, 4000 (top curve). 
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Figure 2(a): 
Legend on graduated horizontal axes: "boundary plane/crystallographic plane 
angle [deg]." 
Legend on graduated vertical axis: "Scholte wave speed [m/s]." 

Figure 2(b): 

Legend on graduated horizontal axes: "boundary plane/crystallographic plane 

angle [deg]." 

Legend on graduated vertical axis: "Scholte wave speed [m/s]." 



Table 1. Values of the elastic stiffnesses (10^° N/m^), density (kg/m^), and 
surface (Rayleigh) wave speed (m/s) for 3 crystals. 



crystal 


Cn 


C22 


Cu 


C16 


C26 


Cm 


Ps 


vr 


ice (-5°C) 


1.38 


1.38 


0.707 








0.3365 


940 


1766 


gypsum 


50.2 


94.5 


28.2 


-7.5 


-11.0 


32.4 


2310 


3011 


terpine 


1.25 


0.99 


0.38 








0.346 


1110 


1644 


germanium 


12.92 


12.92 


4.79 








6.70 


5320 


2936 
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